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We derive an inequality for the distribution of quarks with nonzero orbital angular momentum, and 
thus demonstrate, in a model-independent way, that a nonvanishing anomalous magnetic moment 
requires both a nonzero size of the target as well as the presence of wave function components with 
quark orbital angular momentum L| > 0. 
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O ! I. INTRODUCTION 

(N 

C . Detailed measurements of the spin structure of the nucleon [ESI3;SIEIESI1;0. ^> ^| have revealed that only 
a small fraction of the nucleon spin is carried by the quark spin. This result immediately raised the question, which 
degrees of freedom carry the rest. Unfortunately, both orbital angular momentum and the gluon spin are difficult to 
access experimentally, and therefore little rigorous information exists about quark orbital angular momentum. 
, Meanwhile, many qualitative statements regarding orbital angular momentum have been made. For example, when 
^ ' one expresses the matrix element for the anomalous magnetic moment of the nucleon in terms of light-cone wave 
^\ [ functions (summed over all Fock components), a nonzero anomalous ma gne tic moment can only result when there is a 
nonzero probability that the vector current flips the nucleon helicity |l2lll,'^ . Since the same matrix element conserves 
the spin of the quarks it is evident that some orbital angular momentum must be transfered to the quarks. Hence 
a nonzero anomalous magnetic moment can only occur when the target wave function contains components with 
nonzero orbital angular momentum. While this argument is rigorous, it leaves open quantitative questions regarding 
the norm of those wave function components or perhaps the resulting net L].. Within models, it has also been found 
that a point-like object cannot produce a nonzero anomalous magnetic moment |l2l and within this model one 
can even derive quantitative bounds. 

Similarly, the observation of a nonzero Sivers effect by the HERMES collaboration ^| seems to indicate wave- 
O I function components with nonzero LI since the effect requires an interference between initial nucleon states that have 
1^ opposite helicity. Furthermore, orbital angular momentum seems to play a central role in all models for the Sivers 
function P EI El m in . 

^ 1 The main purpose of this note is to make some of the statements regarding the anomalous magnetic moment more 
quantitative. In order to accomplish this goal, we start from the matrix element that yields the generalized parton 
^ . , distribution and apply the Cauchy-Schwarz inequality which will then provide a lower bound on the norm of wave 
' function components with nonzero orbital angular momentum. 
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II. DECOMPOSITION OF THE NUCLEON SPIN 



In QCD, there is no unique way to decompose the nucleon spin into quark spin, quark angular momentum, gluon 
spin, and gluon orbital angular momentum. For example, Ji has considered a decomposition based on the Af*'^^ 
component of the angular momentum tensor j22l | 

M°^^ ^^q^ (rxDy q+ ^q''a^q + 2Tr£;^ (r x dJ ' A^ + Tr x 1^ . (1) 

The matrix elements of the four terms in Eq. are interpreted as the quark orbital angular momentum, quark 
spin, gluon orbital angular momentum, and gluon spin respectively. One advantage of this decomposition is gauge 
invariance, another is the fact that the matrix elements of the first term in Eq. can be probed in deeply virtual 
Compton scattering. A disadvantage of this decomposition is that the orbital angular momenta in Eq. JQ) contain 
interactions through the gauge covariant derivative. 

An alternative decomposition of the total angular momentum starts from Af in light cone gauge = A'^ +A^ — 
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HI 

= lq\^ (^f X idj ' q+ + ^g^75g+ + 2TyF+^ (r x idj + e+^^^TrF+M^ (2) 

Here g+ = ^7^7^? is the dynamical component of the quark field operators. The obvious disadvantage of using 
Eq. (O to interpret the total angular momentum of the nucleon, which we all know is equal to ^h, as a sum of 
quark and gluon orbital and spin angular momentum, respectively, is the fact that Eq. Q is not gauge invariant. 
However, it is invariant under the residual group of gauge transformations that leave A^ — and where the A^ satisfy 
antisymmetric boundary conditions at x~ = ±oo. The main advantages of using Eq. (0) to decompose the nucleon 
spin is the fact that all terms in Eq. ^ are quadratic in the fields. Unlike Eq. there is therefore no ambiguity 
as to the interpretation of interaction terms. Another advantage is the fact that the matrix elements of the various 
terms in Eq. |(2Jl can be easily expressed in terms of light-cone wave functions. 

We should emphasize that although we have listed here some of the most obvious advantages and disadvantages 
of these two possibilities for decomposing the angular momentum, we do not consider one or the other superior in 
general. However, for specific applications one of these two decompositions may be advantageous. In particular, if 
one wants to place constraints on light-cone wave functions, which find many applications in hadron phenomenology, 
then one may prefer the decomposition of the orbital angular momentum based on M^^"^ |(5J. In the rest of this paper 
we will exclusively study the light-cone decomposition of Jz, based on Eq. 

We thus consider in the following the orbital angular momentum of quarks with flavor q 




and we will investigate its role in nucleon spin-flip matrix elements, such as the anomalous magnetic moment. 

As it stands, Eq. Q is gauge dependent. We therefore need to specify the gauge. First of all, we impose the 
light-cone gauge condition A~^ = 0. However, this does not yet completely fix the gauge since A'^ remains zero under 
gauge transformations with a phase that only depends on xj^ . We thus impose as an additional condition that the _L 
component of the gauge field satisfies anti-periodic boundary conditions at light-cone infinity 

A_L(a;^ = oo,x^) = — A^(a;^ = — oo,xj_). (4) 

With this additional condition, the definition Q of L\ becomes unique. 

Although all results that we will be deriving will be valid in this gauge only, we believe they will nevertheless 
be useful since a lot of hadron phenomenology is based on light-cone wave functions and their use usually implies 
the use of light-cone gauge. It is therefore fair to say that although our results will not have a gauge-independent 
interpretation, they will nevertheless provide additional insights about light-cone wave functions of hadrons. 

III. ANOMALOUS MAGNETIC MOMENT AND GPDS 

We first consider the generalized parton distribution E{x, 0, — A^), which appears in non-forward nucleon spin-flip 
matrix elements of light-cone correlation functions. For purely transverse momentum transfer p+ = p'^ one finds 

(P+A,t|Ona;,0i)|P,T) = i/n.^,0,-Ai) (5) 

(P+A,T|On^,0^)|P,i) = -^^^F'^(x,0,-Ai) (6) 

where 

On.T,bi)= y'^e^f+---g-(0-,bi)7+<7(x-,bx) (7) 

(in gauges other than light-cone gauge = one needs to insert a Wilson line gauge link in Eq. {T))). In the forward 
limit, integration over x yields the anomalous magnetic moment contribution from quarks with flavor q 

Since the operator in Eq. is chirally even, the matrix element is diagonal in quark spin. Since the matrix 
element involves a nucleon spin-flip, angular momentum conservation thus requires a change in quark orbital angular 
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momentum. As a result, the mere fact that the anomalous magnetic moment of the nucleon is nonzero implies that 
there must be components in the nucleon wave function that have a nonzero orbital angular momentum |13| (of course, 
this does not necessarily mean that there is any net orbital angular momentum). In the following we will try to make 
this statement more quantitative, i.e., we will attempt to place constraints on the probability to find nonzero orbital 
angular momentum components in the nucleon wave function. 

In order to facilitate the separation between intrinsic orbital angular momentum and the orbital angular momentum 
due to the motion of the entire nucleon, we first switch to a representation of states that are eigenstates of the transverse 
center of momentum 

|p+,Ri,A) =AfJ d^p^\p+,p^,X)e'''^-'^\ (9) 
where TV is a normalization constant. In this basis, we can define impact parameter dependent parton distributions 



as 



n'^{x,h^) = (p+,Ri,t|0«(x,bx)|p+,R±,T) (10) 

^k-i)'-^ = (p^R.,T|o'(..b,)|p^R,.^>. dD 

The impact parameter dependent PDFs are related to GPDs via a simple Fourier transform [23. l25l 12^ 

n{x,h^) ^ j ^0^e'^--^-H{x,Q,-Al) (12) 



(27r)2 

f(x,bx) = y'^e^^-^-i?(x,0,-Al) 



(13) 



and the normalization is such that / (i'^bj_q(2;, bj_) = q(x). 

Since the GPDs provide simultaneous information about the longitudinal momentum and the transverse position 
of partons, it is instructive to introduce creation operators in this "hybrid space". For example, for x > we define 
destruction operators for quarks with momentum fraction x at transverse position b_L 

b{xp+,h^) = j dx-e'P+''''-'q+{x-,h^) (14) 

and similarly for anti-particles (a; < 0). For now we will suppress the helicity of the quarks and consider only quantities 
that have been summed over the helicities of the quarks. In terms of these, the impact parameter dependent PDFs 
take on a particularly simple form (x > 0) 

n{x,h^) - (p+,R^,T|6t(2;p+,b^)6(a;p+,b^)|p+,R^,T) (15) 

^ f7|- + *7|-V(^'bi) = {p^,^^A\b\xp+,h^)b{xp+,hi_)\p+,B.^,i), (16) 



2M \dh^ dbyy 

which emphasizes their physical interpretation as densities. Upon introducing 

Bl{x)= j (fh^{b^ + iby)b\xp+,h^)h{xp+,h^) (17) 
and integration by parts we find 

Ri, T| B\{x) \p+, Ri, i) = ^i?'(a:, 0, 0). (18) 

We can use this result to provide a formal proof that AL^ = 1 in the matrix element defining E'^{x, 0, 0). Indeed, one 
easily verifies the commutation relation 

[Ll,Blix)]^Blix), (19) 

which proves that has nonvanishing matrix elements only between states that differ by one unit of orbital angular 
momentum L| for flavor q. This observation is consistent with results based on overlap integrals of light-cone wave 
functions |13(|. 
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IV. ANGULAR MOMENTUM DECOMPOSITION OF PDFS 



In order to derive some quantitative constraints on the light-cone wave functions of hadrons, we first introduce an 
angular momentum decomposition for parton distributions, i.e., 

&^(xp+,b^)=^6j„(a;p+,b^) (20) 

m 

where 

hl{xp+,h^) = e+™^^ r *^6t(^p+^bl) (21) 



where hx = |b_L| cos0 and hy = |b^| sincj), and h'^ = \ cos(/)[, and h'y = |b_L| sin^^, respectively. 

Of course, the above decomposition into orbital angular momentum components is scale dependent due to the scale 
dependence of light-cone wave functions |28l |. However, we imagine performing such a decomposition at a fixed scale 
at which we perform the analysis of this decomposition. 

The physical interpretation of the creation and destruction operators H21|l is that they create quarks with m units 
of angular momentum (in the z-direction), i.e., 

[L^, 6j„(xp+, b^)] = mhlXxv^ (22) 

In terms of these angular momentum projected creation operators, the operator B\ appearing in the impact parameter 
space representation of the matrix element for the anomalous magnetic moment takes on the form 

Bl = Y. f '^'^^ {bx+iby)bl^,{xp+,h^)b„Axp+,h^). (23) 

Inserting the angular mode expansion into Eq. H18|l and applying the Cauchy-Schwarz inequality to scalar products 
between states yields 



2M 



E-^ixAO) ^ ^ f d^h^ibx + zby){p+,0^,^\bl^,ixp+,h^)b^{xp+,h^)\p+,0^,l) 
<J2j rf'bi |bi| ^J {p+,0^, Tl bl^,{xp+,h^)b„,+,ixp+, bi) \p+,0^, T) 



{p+,0^,i\bl{xp+,h^)b^{xp+,h^)\p+,0^,i). (24) 



In order to simplify the notation, wc now introduce the distribution of partons with orbital angular momentum m in 
a target with spin f 

ql{x)= y"d2bi(p+,Oi,T|6L(xp+,bi)6„(a;p+,bi)|p+,Oi,r) (25) 
as well as the b^-weighted distribution of partons with orbital angular momentum m 

b^l{x)= J d^h^{p+,0_^,^bl{xp+,h^)b,^{xp+,h^)\p+,0^,^)hl. (26) 

Obviously we have 

qU^) = q^-M (27) 

b'Jix) = b'llix). (28) 
To each term in the sum in Eq. 124|) we now apply the Cauchy-Schwarz inequality for integrals 



J d2biV/(b±).9(b±) <JJ d^h^f{h^)J J d^h^h^). (29) 
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For m > we apply this inequality with 

/(b±) - J d^h^{p+,0^,^\bl^,ixp+,h^)b„,+iixp+,h^)\p+,0^j)=ql^,{x) 

5(b±) = J d^h^hl{p+,0^,l\bl{xp+,h^)b„-,ixp+,h^)\p+,0^,l) =b^J{x), 
while for m < we identify 

f{h±) = J d'h^hl{p+,0^,^bl^,{xp+,h^)b^+,{xp+,h^)\p+,0^,^)^b^j^^,{x) 

9{h±) = J d^h^{p+,0^,l\blixp+,h^)b„,{xp+,h^)\p+,0^,l) ^qLix), 



yielding 



Ei{x,0,0) 
2M 



m<0 



m<0 



r?j.>0 



m>0 



2E V9™+i(^)^-i(^ 



(30) 



(31) 



(32) 



< 2 



\ 



k m>0 



91+1(2;) E^-»(^) 



l,n>0 



where in the last step we used {fm, g-m > 0) 



E V /mSn 



< 



\ 



E/™ Eff" • 

^ m / \ n / 



(33) 



(34) 



If we now introduce the distribution of quarks with positive angular momenta in a nuclcon with spin in the +z 
direction 



(35) 



as well as the contribution to b^(a:) from quarks with zero or negative L| 

6L<o(^) = E^«^(^) 



(36) 



our result can be cast into the form 



fEi{x,0,0) 
V 4M 



< 



qL,>i{x) bl^^oix), 



(37) 



which is the main result of this paper. A slightly stronger inequality can be obtained by keeping track of the quark 
helicity. For this purpose we note that both quark helicities contribute equally in Eq. p8|) . Upon repeating the above 
analysis starting from quarks with helicity |, one thus arrives at 



(38) 



where for example q]^ ^{x) denotes the distribution of quarks with helicity ]. and positive orbital angular momentum 
in a nucleon with helicity f. An even stronger version of Eq. I|38() can be achieved following Ref. 29] by starting with 
7"'"(1 ± 75)72 instead of only 7"*" in Eq. In that case only one of the two terms on the r.h.s. of H38(l remains, and 
gets replaced by |3C| . Additional simplified inequalities can be derived if one neglects quarks with orbital 
angular momenta |L^| > 2. 
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V. 



DISCUSSION 



First of all, our result illustrates that a nonvanishing anomalous magnetic moment (and hence nonvanishing 
i?(x,0, 0)) imphes both 



• a nonvanishing ± size 

The fact that the size of hadrons needs to be nonzero and that nonzero angular momentum components need to be 
present if a state has an anomalous magnetic moment has been observed before in light-cone wave function models of 
hadrons [T^ll4j . What is new in our paper is the fact that our result provides a model independent, quantitative lower 
bound for the _L size distribution, which unlike in Ref. is connected to the quark orbital angular momentum. 
Furthermore, we have shown that there must be quarks with orbital angular momentum in the same direction as 
the nucleon spin — regardless of the sign of the anomalous magnetic moment. This is at first surprising, since one 
might have expected that if the anomalous magnetic moment is due to an orbital angular momentum then, when 
E(x, 0, 0) < the quark should orbit opposite to the nucleon spin. However, one needs to keep in mind that we make 
no statement about the net orbital angular momentum. We only derived a lower bound on contributions to the net 
LI from modes with > 0. 

The necessity of quark orbital angular momentum for an anomalous magnetic moment is particularly surprising 
from the nonrclativistic quark model (NRQM) point of view. In the NRQM the anomalous magnetic moment is 
entirely due to the spins of the quarks and no orbital angular momentum is needed. However, in this point the 
NRQM is not completely consistent: For example, when a massive Dirac particle is in some bound state then, due 
to the localization, the particle must necessarily possess a nonzero momentum and thus a nonvanishing lower Dirac 
component. The lower component for a massive s-wave quark has orbital angular momentum = 1. Since even the 
d quark appears in the NRQM wave functions with both spin up and down, this relativistic effect gives rise to wave 
function components with both positive and negative orbital angular momentum. In our inequality H37|l . the lower 
bound on orbital angular momentum is proportional to the inverse radius squared of the quark distribution. Strictly 
speaking the non-relativistic limit requires Rruq —^ oo and hence RM — > oo. In this limit the lower bound on qL^>i 
goes to zero. Therefore, in a system where the nonrclativistic approximation would be allowed, a nonzero anomalous 
magnetic moment would not necessarily require nonzero orbital angular momentum. 

What is also surprising is that, once E{x, 0, 0) is nonzero, there need to be quarks with positive orbital angular 
momentum — regardless of the sign of E{x,0,0). This is because for E{x,0,0) to be nonzero there ought to be a 
nonzero matrix element between a nucleon with spin up and a nucleon with spin down in Eq. (|18l) . For example, 
if the active quark in the initial state had L| = 0, then it needs to have L| = — 1 in the final state, i.e., L| in the 
same direction as the nucleon spin. Likewise if = in the final state, must have been +1 in the initial state in 
Eq. (|18|l . Similarly one can convince oneself in the general case (initial and final state have nonzero L|) that a nonzero 
E{x, 0, 0) always requires wave function components with > 1 in a nucleon that has Sz =T- This is surprising since 
in Ji's relation the sign of E seems to suggest the sign of L|. However, one needs to keep in mind that Ji's relation 
deals with the net L| and secondly the angular momentum decomposition to which the Ji relation applies ^ does 
not have to be the same as the angular momentum decomposition in the light-cone framework |(5J. 

In order to get some quantitative feeling for the lower bound, we consider a model for w-quarks in the valence region 
where (k — 2) 



where we pick R = 0.5/m for the transverse size, i.e., MR = 2.5 and we let b'^{x) cx (1 — x)'^q{x), consistent with a 
finite size for large x. Since no data is available on E{x, 0, 0) and 6^(x), we simply make an educated guess regarding 
these functions in order to be able to make an order of magnitude estimate for the probability to find nonzero orbital 
angular momentum. The ansatz for E(x, 0, 0) is motivated by the constraint that E(x, 0, 0) must vanish by two powers 
of (1 — x) faster than q(x) for large x. The ansatz for b^{x) was motivated by the constraint that the nucleon has 
a finite transverse size for large x. The rest of the above ansatz was guided by the ad hoc requirement to make the 
model as simple as possible and therefore the idea is that this model should merely serve as a guide for what order of 
magnitude one should expect ( "back of an envelope estimate" ) with parameters that are within the range of what is 
expected in the valence region. With the above model parameters we find for the distribution of quarks with positive 
orbital angular momentum 



• wave function components with positive quark orbital angular momentum (for Sz = +^) 



q{x) 
E{x,0,0) 
M^b^ix) 



8(1-2;)=^ 

6k{1-x)^ = 12(1 -x)^ 
M'^R^8{1 - x)^ 



(39) 
(40) 
(41) 



^{l-xf<qL.>i{x), 



(42) 
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i.e., only some % compared to the distribution summed over all orbital angular momenta. 

While this is much smaller than current estimates based on Ji's relation (see for example Ref. [sj] for an up to date 
estimate) , one should keep in mind that the orbital angular momentum obtained from the Ji relation differs from the 
orbital angular momentum from the light-cone decomposition (j^J. 

Although our lower bound is not very spectacular, it still provides the first model-independent lower bound on the 
distribution of quarks with positive angular momentum distribution. One may also wonder how such a low bound 
is consistent with transverse flavor dipole moments w 0.2/m for a moving nucleon that are quite significant and 
almost the same order of magnitude as the _L size of the nucleon. The important point here is that the _L distortion, 
which is described by our starting equation E, contains a piece involving the overlap between the L| = component 
of the wave function and the = 1 component, i.e., it is linear in the amplitude for finding a quark with nonzero 
orbital angular momentum, whereas the probability for L| > is quadratic in that amplitude. The situation is thus 
somewhat analogous to the quadrupole moment of the deuteron, which is linear in the d-wave component and quite 
large, while the actual d-wave probability is tiny. 

VI. SUMMARY 

We have derived an inequality, which provides a model-independent lower bound on the norm of wave-function 
components with nonzero quark orbital angular momentum. Although the numerical bound thus obtained is not 
very strong — requiring only a few % probability for nonzero orbital angular momentum - our result represents the 
first quantitative estimate for wave- function components with nonzero L|. Moreover, the resulting constraints on the 
light-cone wave functions of the nucleon, which only enter the distributions quadratically, arc much stronger. The 
bound that we derive involves the transverse size (rms-radius) times the mass of the nucleon. This quantifies the 
known result that a point-like particle cannot have a nonzero anomalous magnetic moment. Unlike L], itself, which 
depends quadratically on wave- function components with nonzero L^, the GPD _B(a;, 0,0) requires only that is 
nonzero on one side of the matrix element. This illustrates why we were only able to derive a bound on the probability 
for nonzero that depends quadratically on E{x, 0, 0), which is why our bound is so weak. In contradistinction 
obtained through the Ji relation would involve E{x, 0, 0) linearly. Nevertheless, we hope that the constraints derived 
in this work will be of use in developing better models for light-cone wave functions describing these interesting 
observables. 
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